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ABSTRACT 

Studying the reflection of X-rays off the inner edge of the accretion disc in a neutron star low- 
mass X-ray binary, allows us to investigate the accretion geometry and to constrain the radius 
of the neutron star. We report on a NuSTAR observation of 4U 1608-52 obtained during a faint 
outburst in 2014 when the neutron star, which has a known spin frequency of u = 620 Hz, 
was accreting at ~l-2 per cent of the Eddington limit. The 3-79 keV continuum emission 
was dominated by a T ~ 2 power law, with an ~l-2 per cent contribution from a fcTbb—0.3- 
0.6 keV blackbody component. The high-quality NuSTAR spectrum reveals the hallmarks of 
disc reflection; a broad iron line peaking near 7 keV and a Compton back-scattering hump 
around ~20-30 keV. Modelling the disc reflection spectrum points to a binary inclination of 
i ~30°-40° and a small ‘coronal’ height of h <8.5 GM/(?. Furthermore, our spectral anal¬ 
ysis suggests that the inner disc radius extended to i?in—7-10 GM/c^, close to the innermost 
stable circular orbit. This constrains the neutron star radius to R <21 km and the redshift from 
the stellar surface to z >0.12, for a mass of M = 1.5 Mq and a spin parameter of a = 0.29. 

Key words: accretion, accretion discs - stars: individual (4U 1608-52) - stars: neutron - 
X-rays: binaries. 


1 INTRODUCTION 

Measuring the mass and radius of neutron stars would constrain 
their equation of state and thereby give insight into the behaviour 
of matter at supranuclear densities, a challenge that cannot be 
achieved in terrestrial laboratories. Since many proposed equations 
of state predic t a range of masses for a given radius (see e.g., 
iLattimeijboi ll . for a recent review), it is of particular interest to 
measure neutron star radii. Low-mass X-ray binaries (LMXBs) are 
promising targets to obtain such information, as the hot surface of 
the neutron star may be directly visible in the X-ray band. 

When residing in an LMXB, a neutron star accretes mat¬ 
ter from an accretion disc that is fed by a low-mass (< 1 M©) 
companion star. Often LMXBs are transient, exhibiting orders 
of magnitude variation in their X-ray luminosity driven by simi¬ 
larly large variations in the mass-accretion rate. During quiescent 
episodes the X-ray luminosity is low (Lx — erg s“^) 

and little matter is thought to be accreted on to the neutron star. 
Thermal emission from the stellar surface may then be detected 
with sensi tive X-ray instrumen t s and used to measure neutron star 


radii (e.g..[Rutledge et alJ 19991: 

Webb & Barretl2007 

:lGuillot et al.l 

|201 lLl2013l;IServillat et al.ll2012 

; Heinke et all 20061 

2014). 
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During outburst phases, on the other hand, matter is rapidly 
accreted on to the neutron star, typically generating an X-ray lu¬ 
minosity of 1-100 per cent of the Edd ington limit (Lsdd — 
3.8 X 10^® erg s~^: fkuulkers et alj|200^ . Although the overall 
X-ray emission is then dominated by that of the accretion flow, 
the neutron star becomes visible during type-I X-ray bursts; bright 
flashes of X-ray emission resulting from unstable thermonuclear 
burning of accreted H/He on the surface of the neutron star. Mod¬ 
elling the black-body s pectra of these events also facilitates ra- 


dius measurements (e.g. .Ivan Paradiisll 19791: 

Fuiimoto & Gottwald 

1989 

: iGiiver et alj 120101: Suleimanov et alj 

201 1|; 

Poutanen et al. 


2014 . 


Both these methods rely on using appropriate neutron star 
atmosphere models, which has spurred in tense discussion (e.g., 
ISuleimanov et alJl201ll : [Heinke et al.ll20l4 . Alternative means of 
constraining neutron star radii are offered by studying X-rays re¬ 
flected off the inner edge of the accretion disc, which may ex¬ 
tend very close to the stellar surface (e.g., ICackett et al.l 1201^ ; 
iMiller et aDl2013^ . The observable effects of disc reflection are a 
broad emission line in the Fe-K band (6.4-6.97 keV) and a Comp¬ 
ton b ack-scattering hump peaking at ~20-40 keV l lFabian et al.l 
Il989l) . 

Studying reflection spectra also offers valuable insight into 
the accretion geometry, such as the inner radius and inclination 
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Figure 1. NuSTAR FMPA/FMPB summed, background-corrected light 
curve at 1-s resolution (3-79 keV). The time of the SwiftfXKl observation 
is indicated by the horizontal bar. 


of the accretion disc as well as the height of the illuminating 
X-ray source, and how this is affected by the accretion rate or 
the magnetic field of the neutron star. Analysis of broad Fe lines 
in several neutron star L MXBs has revealed inner disc radii of 
J?in— 5-20 GM/c^ (e.g., CacketteLal l201(]l : lEgron et all 1201 ll : 
Is anna et al ]|2014l :lDi Salvo etal.1120151 for recent studies), i.e.. 


close to the innermost stable orbit (ISCO) in the Schwarzschild 
metric (7?isco—6 GM/c^). However, in some LMXBs the disc 
appears to be truncated at larger radii due to the pressure exerted 
by the magnetic fi eld of the neutron star (e.g., GRO J1744-28 with 
J?in—85 GM /: lOegenaM et alJ|20l4 , or due to evaporation of 
the inner disc at lo w accretion rates (e .g., HETE J 1900. 1-2455 with 
7?in~25 GM/c^ : fp^itto et al.ll2013t) . 

Detecting a Compton hump in addition to a broad iron line 
requires high sensitivity at energies > 1 0 keV, such as provi ded by 
the recently launched NuSTAR satellite jHarrison et alj2013l) . High 
quality NuSTAR spectra ca n allow for new vi ews of the accretion 
geometry in LMXBs (e.g., iMiller^t^ 2015ll. and constraints of 
neutron star radii (e.g., Ser X-h lMiller et alJl2013al) . 


2 OBSERVATIONS AND ANALYSIS 

We observed 4U 1608-52 simultaneously with NuSTAR and Swift 
JOehrels et al.ll200^ on 2014 October 16-17. Data reduction and 
analysis was carried out using tools incorporated in HEASOFT ver. 
16.6. Throughout this work we assume a distance of D=3.6 kpc 
and report errors as 90 per cent confidence levels. 


2.1 NuSTAR 

NuSTAR observed 4U 1608-52 between 23:00 UT on 2014 Octo¬ 
ber 16 and 17:10 UT on October 17 (ID 90002002002). Standard 
screening and processing with NUSTARDAS (ver. 1.4.1) resulted in 
~32 ks on-target exposure time. We created light curves and spec¬ 
tra for the FPMA and FPMB employing the NUPRODUCTS tool. To 
obtain source events we used a circular extraction region with a ra¬ 
dius of 120 arcsec, whereas a region of the same dimensions placed 
away from the source was used for the background. 

Light curves obtained for the two modules were first back¬ 
ground subtracted and then summed using LCMATH. Initial inspec¬ 
tion of the spectra revealed that the source was detected signif¬ 
icantly above the background in the entire NuSTAR band (3-79 
keV), and that the separate FPMA/FPMB data showed excellent 
agreement. We therefore combined these using ADDASCASPEC, 
which also generates the combined background spectrum and ancil¬ 
lary response file (arf). A combined redistribution matrix file (rmf) 
was created via ADDRMF, weighting the responses of the two mod¬ 
ules by their exposure times. Using GRPPHA we grouped the data 
to a minimum of 20 photons per spectral bin. 


2.2 Swift 

Simultaneous Swi/i/XRT data were obtained to provide energy cov¬ 
erage down to ~0.5 keV. 4U 1608-52 was observed for ~L7 ks 
between 23:56 UT 2014 October 16 and 01:40 UT October 17 
(ID 32322019), with the XRT operated in windowed timing (WT) 
mode. Using XSELECT we extracted source events from a box of 
120 arcsec long and 40 arcsec wide. A box of the same dimen¬ 
sions placed away from the source was used to extract a background 
spectrum. 

The WT light curve showed a stable intensity along the obser¬ 
vation (~5 counts s“^), and we therefore extracted a single aver¬ 
age spectrum. An arf was created using XRTMKARF and the latest 
rmf (ver. 15) was sourced from the calibration data base. The spec¬ 
tral data were grouped to a minimum of 20 photons per bin. 


I.I 4U1608-52 


In this work we report on the disc reflection spectrum measured 
by N uSTAR for the neutron star 4U 1608-52 jGrindlav & Gursl^ 
ll976l : lThnanbaum et al.lll97^ . This transient LMXB is frequently 
active, with a ccretion outbursts typically recurring once every 
c::l-2 yr (e.g ., Lochner & Roussel-DuDrellT994 IChen et al.ll 19^ : 
ISimonll2()04 ; I Galloway et alj 2008ll . Tyre-1 X-ray bursts are regu¬ 
larly observed and h ave allowed for a dis t ance estimate of D ~ 
2.9 — 4.5 kpc (e.g.. iGallowav et ^l2008l : [Poutanen et al.ll20l3l . 
Moreover, rapid oscillations detected during type-I X-ray bursts 
revea l ed that the neutron s tar spins aX v = 620 Hz l lMuno et al.l 
1 200 ll : [galloway et al.ll200^ . Renewed ac tivity was detected from 
the source on 2014 October 5 with MAXI jNegoro et al.ll2014h . 


3 RESULTS 

3.1 NuSTAR light curve 

Fig.[T]shows the light curve of 4U 1608-52 obtained with NuSTAR. 
The source was detected at an average intensity of ~70 counts s“^ 
(3-79 keV, two modules combined). A strong ~200-s long increase 
in intensity was registered ~12 ks into the observation, caused by 
the occurrence of a type-I X-ray burst from the source. Analysis of 
this event will be presented in a separate work. 


3.2 Spectral continuum 

We examined the spectral continuum by fitting the NuSTAR data 
(excluding the X-ray burst) together with the 5wi/f/XRT data in 
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Figure 2. Unfolded NuSTAR (black) and Swift (red) spectra. The solid lines 
indicate fits to an absorbed, phenomenological continuum consisting of a 
r ~ 2.0 power law (dashed lines), a cool ~ 0.1 keV blackbody (dot¬ 
ted curve), and a hotter /cTbb — 0.4 keV blackbody (dash-dotted curves). 
The bottom panel shows the data-to-model ratio. 
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Figure 3. NuSTAR data-to-model ratio in the Fe-K line region for a spectral 
continuum without disc reflection (rebinned for visual claiity). 


XSPEC (ver. 12.8: lArnaudlll996h . A constant multiplication factor 
was included to account for calibration differences. To model the 
interstellar absorpt ion we used TBABS wit h VERN cross-sec tions 
jVerner et alJll996h and WILM abundances iWilms et al.ll2Q0(]|) . 

Any realistic combination of continuum models left large pos¬ 
itive residuals around :^5-8 and 20-30 keV. This is illustrated by 
Fig.m where we show a fit to a continuum consisting of a power 
law and two blackbody components. The prominent residuals can 
be interpreted as a broad Fe-K emission line, shown in more de¬ 
tail in Fig.[^ and the corresponding Compton hump. We therefore 
proceeded by modelling our data with physical reflection models. 


3.3 Reflection spectrum 

As can be seen in Fig. the continuum emission at energies of 
>2 keV is strongly dominated by a hard spectral component. We 
therefore chose the RELXILLLP model, which calculates disc reflec- 


tion features due to an illuminating power-law source Joi 


arcia et 


aO 


|2014) . This new model features higher spectral resolution and up¬ 
dated atomic data with regard to other reflection models. More- 
over, it combines th e relativistic convolution kernel RELCONV_LP 
I Dauser et al .I 2 OIOI) with the reflection grid XILLVER dOarcfa et al.l 
l2013h . in which the reflection spectrum is calculated for each emis¬ 
sion angle rather than averaged. Both the illuminating power law 
and the reflected emission are fitted self-consistently. 

In RELXILLLP the hard X-ray source, usually referred to as the 
corona, is assumed to be a point source located at a height h above 
the accretion disc plane along the spin axis of the compact object. 
Although this ‘lamppost geometry’ is an oversimplification (albeit 
compact, the corona is likely not a point sou rce), it provides an 
adequ ate description for several LMXBs (e.g.. lMiller et alJboiSbl 
l2015h. The emissivity profile is calculated based on the height h 
rather than assuming a (broken) power-law profile. 

The fit parameters of the RELXILLLP model are the dimen¬ 
sionless spin a, the binary inclination i, the inner and outer disc 
radii Rin and Rout (expressed in terms of i?isco for the given a), 
the ionization parameter log 5, the iron abundance Ape (with re¬ 
spect to Solar), the reflection fraction i?refl (ratio of the reflected 
to primary emission between 20 and 40 keV), the normalization 
Arefl, and the index F and high-energy cutoff Ecut of the power 
law. Since 4U 1608-52 has a known spin frequency, the dimen¬ 
sionless spin parameter can be calculate d. For neutron st ars this 
can be approximat ed as a ~ 0.47/Prms1 ^ Braje_et_ay2000|). Given 
that u = 620 Hz jMuno et al.ll200lHGallowav et alj|2008l) . we set 
a = 0.29 in all our fits. Moreover, we fixed Rout ~ 400 F?isco 
since the emissivity profile drops off steeply with increasing radius, 
so that reflection fits are not sensitive to the outer disc radius. 

Fitting the spectral data required one or two soft emission 
components in addition to the reflection model, for which we used 
the model BBODYRAD. The thermal emission may originate from 
the neutron star surface, the accretion disc, or a boundary layer be¬ 
tween the two. We found that the relatively high extinction towards 
4U 1608-52 and the limited number of counts in the Swift data 
caused considerable uncertainty in the spectral shape at <3 keV. To 
investigate to which extent this influenced the (reflection) spectrum 
at higher energies, we explored a number of different fits. These are 
summarized in Table[T]and discussed in more detail below. 

We initially fitted the combined NuSTAR and Swift data using 
the full 0.5-79 keV energy range, which required two soft compo¬ 
nents (model 1 in TablelTJ. However, the obtained normalization of 
the ~0.1 keV blackbody would imply an unphysically large emis¬ 
sion radius of >500 GM/(?. Possibly, the requirement for this soft 
component is due to a calibration uncertainty that may arise for 
highly absorbed sources observed in WT mode and causes a bump 
at energies of ~0.5-l kevQTo mitigate these possible effects we 
also performed fits in the 1-79 keV range, which required only one 
blackbody component (model 2 in Table[T). 

We also explored fits using only the NuSTAR spectral data. 
When leaving the hydrogen column density free (model 3 in Ta¬ 
ble [T]l the obtained value is notably larger than for fits 1-2, and 
also larger than the value s of N-g ~ ( 0.5 — 2) x 10^^ cm“^ 
reported previous ly (e.g., IPenninx et alJll98^ : iKeek et Tl 1200Sl : 
iGtiver et al. I l20i^ . We therefore also fitted the NuSTAR spectrum 
with A^h = 2 X 10^^ cm“^ fixed (model 4 in TableO, shown in 
Fig.lH 

Importantly, we find that despite the considerably different pa- 


^ http://www.swift.ac.uk/analysis/xrt/digest_cal.php 
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20-30 keV. We fitted the spectral data with the versatile disc reflec¬ 
tion model RELXILLLP to investigate the accretion geometry of 4U 
1608-52, and to obtain constraints on the radius of the neutron star. 


4.1 Accretion geometry of 4U 1608-52 

It is generally thought that in soft X-ray spectral states, when the 
accretion rate is high (>10 per cent of the Eddington limit), the ac¬ 
cretion disc extends to/near the ISCO, whereas in quiescence the 
inne r disc radius lies far from the compact object (>100 GM/c?\ 
e.g., lEsin et al ] [T^ . Therefore it is expected that somewhere 
in between these extremes, in the hard X-ray spectral state, the 
inner edge of the disc should retreat. However, there is consid¬ 
erable debate whether the accretion disc is truncated through¬ 
out hard X-ray spectral states, or whether its inner radius only 
starts to recede below a certain accretion rate, because radii in- 


Figure 4. Unfolded NuSTAR spectrum and model fit (number 4 in Tabled 

opposing results (e.g.. Rvkoff et alJ 2007; Gierlinski et al.l 

2008; 

solid line) consisting of an illuminating power law and resulting disc reflec- 

Tomsick et alj 20091; Reis et al. 201C 

; Done & Diaz Trigd 

201C; 

tion (RELXILLLP; dashed line) plus a blackbody (BBODYRAD; contributing 

Kolehmainen et alj |2014|; IPlant et alj 

20141). This discussion has 


only at <5 keV). The bottom panel shows the data-to-model ratio. 


rameter values that we obtain for the absorption and thermal emis¬ 
sion component(s), the results obtained for the hard X-ray contin¬ 
uum and the reflection spectrum are robust (Table[TJ. Our fits point 
to a moderately low inclination of i~30-40°, consistent with ob¬ 
servations of the optical counterpart (see Section IdTt . The obtained 
inner disc radius lies near the ISCO, i?in—l.3-2.0 i?isco, and the 
coronal height is small, h <8.5 GM/c?. We find a reflection frac¬ 
tion of iirefl —1 and a high ionization parameter of log ^ ~ 3.9, as 
may be expected given the breadth of the Ee line (Fig.j^. 

We note that our spectral fits favour a high value for the power- 
law cutoff; Ecut > 300 keV. Although this is not unprecedented for 
neutron stars, it is well outside the NuSTAR bandpass. This could 
perhaps indicate that the spectral shape of the illuminating X-ray 
source (slightly) deviates from the cutoff power law that is incor¬ 
porated in the RELXILLLP model. 

Our different fits yield an unabsorbed flux in the 3-79 keV 
band of Fz-iq ~2.0 x 10“® erg cm“^ s“^, of which ~l -2 per 
cent is due to thermal emission. The uncertainty in the low-energy 
spectrum causes some spread in the inferred 0.5-79 keV fluxes; 
we obtain To. 5-79 ~(3.0 — 5.2) x 10~® erg cm“^ for fits 
2-A, and ~1.1 x 10“® erg cm“^ for fit 1. The latter may be 
overestimated due to a possible calibration issue in the WT data. 
We therefore estimate that 4U 1608-52 accreted at ~l-2 per cent 
of the Eddington limit during our observations. 


4 DISCUSSION 

The frequently active transient neutron star LMXB 4U 1608-52 
started a new, faint outburst in 2014 October. We obtained simul¬ 
taneous NuSTAR and Swift observations during its initial hard X- 
ray spectral state. The continuum emission in the 3-79 keV NuS¬ 
TAR band was dominated by non-thermal emission that can be de¬ 
scribed by a r~2 power law, with a small (~l-2 per cent) contri¬ 
bution from a soft emission component that can be modelled with 
an ~0.3-0.6 keV blackbody. The hard X-ray continuum and the 
excellent sensitivity of NuSTAR at energies >10 keV resulted in 
a high-quality disc reflection spectrum, showing both a broad Fe- 
K emission line peaking near 7 keV and a Compton hump around 


mainly focused on black hole LMXBs. The accretion geometry 
could possibly be different in the hard X-ray spectral states of neu¬ 
tron star LMXBs, as the stellar surface and anchored magnetic field 
may come into play. 

We estimate that 4U 1608-52 accreted at ~ 1-2 per cent of the 
Eddington limit during our observations, which thus probe a rel¬ 
atively low accretion regime for disc reflection studies of neutron 
star LMXBs. Our spectral fits consistently point to an inner disc 
radius that lies close to the ISCO: i?in=L3-2.0 i?isco- Given that 
77isco—5.05 GMjS for a neutron star spinning at o=0.29, this 
would correspond to R\^= 7-10 G'M/c^ = 15-21 km for a mass of 
M=1.5 M( 7 ) (a reasonab le choice based on the recent overview of 
iLattimer & Steinerll2014l) . Our obtained inner disc radius for 4U 
1608-52 is within the range foun d for several other neutron star 
LMXBs (i7in—5-20 GM/c^; e.e.. lCackett et al.l2010l ; l^ron et al.l 
l201lh . Our study does not indicate significant trunc ation of the ac¬ 
cretio n disc, despite the low accretion rate (see also lPi Salvo et all 
l2015h . Due to uncertainties in our data at <3 keV, we cannot obtain 
reliable radius estimates from the thermal emission in 4U 1608-52. 

Our spectral analysis further points to a moderately low disc 
inclination angle of i ~30°^0°. This is consistent with the lack of 
dips/eclipses from 4U 1608-52 an d with the possible detection of a 
‘superhump’ dWachter et al.ll2002h . which should be best observed 


at low inclination (e.g., Haswell et alJl200lh . Furthermore, we find 


that the illuminating hard X-ray source is likely located close to the 
neutron star, at a height of /i<8.5 GM/c?. This may be consistent 
with the growing consensus that the hard X-ray corona in similar 
accreting systems ( black hole LMXBs and active galactic nu clei) is 
very compact (e.g., iReis & Milleill2013l;lFabian et alj|20l4 for re¬ 
cent discussion). Alternatively, the very small height inferred from 
our reflection fits could point to the boundary layer between the 
accretion disc and the stellar surface as the primary so urce of the 
illuminating hard X-rays (e.g.. lGierlinski & Donell20o3) . 


4.2 Neutron star radius constraints 

Since the disc must truncate at the surface of the neutron star 
if not at larger radii, reflection modelling c an be used to place 
constraints on the neutron star radius (e.g., ICackett et al.ll 201 ^ ; 
iMiller et"^ l2013ah . For a gravitational redshift of 1 -f z = 
l/i^l — 2GM/Rc^, the inner disc radius implied by our fits 
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Table 1. Results from modelling the spectral data with disc reflection. 


Model number 

1 

2 

3 

4 

Fitted data 

Swift+NuSTAR 0.5-79 keV 

Swift+NuSTAR 1-79 keV 

NuSTAR 3-79 keV 

NuSTAR 3-79 keV 

CONSTANT*TBABS’* 

(BBODYRAD+BBODYRAD+RELXILLLP) 

(BBODYRAD+RELXILLLP) 

(BBODYRAD+RELXILLLP) 

(BBODYRAD+RELXILLLP) 

C . 

1.04 ±0.02 

1.02 ± 0.01 

- 

- 

Nil (Xl 022 cm- 2 ) 

2.6 ± 0.2 

1.7±0.1 

3.9 ± 1.1 

2.0 fix 

fcTbbi (keV) . 

0.46 ±0.02 

0.56 ±0.03 

0.43 ± 0.04 

0.31 ±0.03 

Nbbi (km /10 kpc)^ 

6 . 0 I 2 2 X 10= 

l-dlo 4 X 10 ^ 


1.i1q7 X 10'* 

fcTbbz (keV) . 

0.13 ±0.01 

- 

- 

- 

Afbb 2 (km /10 kpc )2 

9 . 9+5 g X 10® 

- 

- 

- 

h(GM/A) . 

1 q+ 6-6 

-*-•^-1.9 

1 q+3'3 
-*-•^-1.9 

0 . 8 t®j 

-*-•^-1.9 

*(°) . 

36.1 ± 6.9 

31.1 ± 6.5 

41.1 ±6.1 

38.8 ± 6.0 

Riu (Xiiisco) • ■ ■ • 


T 7 + 0.6 
-*-• '-0.7 

1 q+0.8 

^•'^-1.3 

1 Q+1-9 

r. 

1.96 ± 0.02 

1.97 ± 0.01 

1.97 ± 0.01 

1.95 ± 0.01 

log?. 

4.0 ±0.2 

4.1±0.1 

3.9 ± 0.2 

3.9 ± 0.1 

Ape (x Solar). 

2.2 ±0.5 

2.5 ±0.7 

1.7 ± 0.7 

1.8 ± 0.4 

Bout (keV). 

4391^?^ 

r;54+l00 

O0+_i5o 

K-| 0+254 

416t“® 

-^refl . 

1.3±0.7 

1.5 ± 0.5 

1.0 ±0.4 

1.0 ±0.3 

Nrefl . 

0.15 ±0.03 

0.14 ±0.03 

0.18 ±0.03 

0.18 ±0.03 

xl (dof) . 

1.18 (1847) 

1.22(1811) 

1.06(1280) 

1.07 (1281) 


Notes. - The constant C was fixed to 1 for the NuSTAR data and left free for the Swift spectrum. The outer disc radius of the RELXILLLP spectral component 
was always fixed at Rout = 400 Risco 2 nd the spin parameter was set to a = 0.29 to correspond to the known spin frequency of u = 620 Hz. The 
reflection fraction Rrefl is calculated between 20 and 40 keV. Quoted errors reflect 90 per cent confldence levels. 


would constrain the neutron star radius to R<21 km, hence the 
gravitational redshift to z:>0.12 for an assumed mass of M = 
1.5 Mq. These constraints from the disc reflection spectrum of 4U 
1608- 52 are consistent with those obtained fr om its type-I X-ray 
bursts jGiiver et alJl2010l : IPoutanen et alj2014h . 
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